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In this work we are reporting on the measurement of the proton-air inelastic cross section 
using the Telescope Array (TA) detector. Based on the measurement of the the proton-proton 

cross section cTp-p value is also determined at = 95lg TeV. Detecting cosmic ray events at ultra 
high energies with Telescope Array enables us to study this fundamental parameter that we are 
otherwise unable to access with particle accelerators. The data used in this report is the hybrid 
events observed by the Middle Drum fluorescence detector together with the surface array detector 
collected over five years. The value of the is found to be equal to 567.0 ± 70.5[Stat.]lt25 [Sys.] 

mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and 
Block, Halzen and Stanev QCD inspired fit and is found to be equal to I 70 I 4 I [Stat.]lJ;®[Sys.]mb. 


I. INTRODUCTION 

Measuring the proton-air inelastic cross section 
from cosmic rays at ultra high energies allows us to 
achieve knowledge of a fundamental particle property 
that we are unable to attain with measurements at cur¬ 
rent accelerators. At the time of writing of this paper the 
highest proton-proton center of mass energy that could 
be attained by the modern accelerators is ^ 14 TeV by 
the Large Hadron Collider (LHC). However, Ultra High 
Energy Cosmic Ray (UHECR) experiments have been re¬ 
porting on the proton-air inelastic cross section starting 
with the Fly’s Eye in 1984 at = 30 TeV [I| and ending 
with the most recent result of the Auger experiment at 
v/s = 57TeV in 2012 0. 

The current high energy models agree in their pre¬ 
dictions of rising proton-air cross section with energy. 
The high energy models are in reasonable agreement at 
lower energies, below 10^® eV, where they are tuned to 
measurements of multi-particle production provided by 
particle accelerators. However the high energy models 
diverge in describing fundamental parameters such as 
hadronic cross sections, elasticity, and secondary parti¬ 
cle multiplicity above 1 PeV where the models rely solely 
on theoretical expectations 0. Studying the energy de¬ 
pendence of the proton-air cross section is important in 
constraining the extrapolation of the hadronic models to 
high energy. 

Detecting UHECR showers provides the opportunity 
to study fundamental particle properties. Optimally, to 
measure the directly, we observe the first point of 

the proton-air interaction slant depth Vi and fit the dis¬ 
tribution of Xi to recover the interaction length Ap_air- 
However, since the observation of the first point of inter¬ 
action to obtain the nucleon-air cross section is not fea¬ 
sible, the inelastic proton-air cross section is calculated 
using the distribution of the observed shower maximum 
Xmax- Ap_air and Consequently cr^bj., are derived from 
the Xmax distribution’s exponential tail. 

In this work, we report on the measurement of the 
proton-air inelastic cross section tTp^b^ using the Tele¬ 
scope Array’s Middle Drum detector together with the 
Surface array Detector (MD-SD) in hybrid mode data 0. 
The method used in this calculation is “the if-Factor 
method”, were the underlying assumption is a propor¬ 
tionality between the tail of the Xmax distribution and 
Xi. Details of the method, the result, and the system- 
atics of the measurement are presented in this work. In 


addition, the Telescope Array proton-air inelastic cross 
section is compared to previous results. Furthermore, 
the proton-proton cross section (Tp_p is calculated us¬ 
ing Glauber theory together with the BlockMdalzen, and 
Stanev (BHS) QCD-inspired fit ( 0 , 0 , 0 , 0 ). The 
proton-proton cross section is also compared to previous 
tTp_p experimental results. Finally, we discuss the sum¬ 
mary and the outlook. 


II. DATA TRIGGER, RECONSTRUCTION, 
AND SELECTION. 

The data used in this analysis is collected by the Tele¬ 
scope Array (TA) detector located in the southwestern 
desert of the State of Utah. TA is an UHECR detec¬ 
tor composed of three Fluorescence Detector (ED) sites 
and the Surface Detector array (SD) 0 as shown in Fig¬ 
ure [H The SD array occupying 700 km^ and is bounded 
by the FDs. The northernmost Fluorescence detector is 
referred to as Middle Drum (MD), while the other two 
southern FDs are named Black Rock Mesa (BRM) and 
Long Ridge (LR) [ig. Moreover, a Central Laser Facil¬ 
ity (CLF) to monitor the atmosphere and calibrate the 
detector is deployed in the middle of the detector and is 
located equidistant from the three FDs. 

The two southernmost detectors, LR and BRM, consist 
of 12 telescopes each, while the MD detector separated 
from the northern edge of the SD by 10 km consists of 14 
telescopes each of which uses a 5.1 m^ spherical mirror. 
The fluorescence light from each mirror is collected to a 
camera containing 256 Photo Multiplier Tubes (PMTs) 
tightly spaced. Seven of the 14 mirrors view 3°- 17° 
in elevation while the rest view 17°- 31°, with a total 
azimuth of 112° between southwest and southeast. On 
the other hand, the SD is composed of 507 scintillation 
counters each 3 m^ in area. The SD scintillation counters 
are spaced on a 1.2 km grid. 

The data used in this analysis consists of MD-SD hy¬ 
brid events. The MD and SD trigger independently. Of¬ 
fline, a hybrid data set is formed by time-matching the 
events from the two detectors. In the monocular mode, 
an event trigger is recorded by the SD when three adja¬ 
cent SDs observe a signal greater than 3 Minimum Ioniz¬ 
ing Particles (MIPs) within 8 /iS. When a trigger occurs 
the signals from all the SDs within ± 32 /xs with ampli¬ 
tude greater than 0.3 MIP are also recorded. Moreover, 
a telescope event trigger for the MD is recorded when 
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two subclusters of the 256 PMT cluster triggered within 
25 ns. Here a subcluster is defined as a (4 x 4) 16 PMTs 
within the 256 PMT cluster. Each subcluster reports a 
trigger when three tubes in that subcluster trigger within 
25 ^s, two of which are adjacent. Finally, multiple tele¬ 
scope event triggers within 100 /rs would be combined 
into a single MD event. 

Events detected by both detectors (MD and SD) within 
2 ^s are combined into one hybrid event. The com¬ 
bined data set with MD and SD time matched events 
are then reprocessed using information from both detec¬ 
tors. The ED overlooking the sky above the SD array 
provides the longitudinal profile of the shower. Mean¬ 
while, the SD provides the event shower core, particle 
density, and hence improves the geometrical reconstruc¬ 
tion significantly. Reference [3] describes the detector 
monocular and hybrid reconstructions of the triggered 
events in more detail. 

To achieve the best X^ax resolution, a pattern recog¬ 
nition technique was applied which selected events with 
a well-defined peak in the fluorescence light profile. This 
technique is described more completely in Reference [^. 
Briefly, each shower profile’s shape was approximated by 
a set of right triangles, and a set of cuts on the properties 
of these triangles was used to reject X^ax events with a 
“flat” profile or an indistinct peak. As shown in Refer¬ 
ence d I, data to Monte Carlo comparison studies showed 
good agreement in basic air shower distributions such as 
zenith angle, azimuthal angle, and impact parameter af¬ 
ter these pattern recognition cuts were applied. 

The data used in this analysis is the MD-SD hybrid 
events collected between May-2008 and May-2013. After 
applying the pattern recognition cuts to the data we are 
left with 439 events. The energy range for this data set is 
between 10^®'^ and 10^®'® eV. With an average energy of 
1018.68 gY^ jg equivalent to a center of mass energy of 
= 95 TeV. Finally, the Xmax resolution of this data 
set achieved after applying the pattern recognition cuts 
is ^ 23 g/cm^ i. 


III. ANALYSIS 


In this paper we determine the value of using 

the AT-Factor method. This method infers the attenua¬ 
tion length and hence the cross section value from the 
exponential tail of the Xmax distribution. This is assum¬ 
ing that the tail of the Xmax distribution is comprised 
of the most penetrating/lighter particles (protons). The 
tail of the Xmax distribution is fit to the exponential 
exp( )^ where A^ is the attenuation length. Km is 

proportional to the interaction length Ap_air: 


Km — ATAp—air — AT 


14.45mp 

p—air 


( 1 ) 


where AT is dependent on the shower evolution model. 
The departure of AT from unity depends on the pion 


TAMAP 



SD East/West Column Number [1200m units] 


FIG. 1. The Telescope Array detector configuration. The 
filled squares are the 507 SD scintillators on a 1.2 km grid. 
The SD scintillators are enclosed by three fluorescent detec¬ 
tors shown in filled triangles together with their field of view in 
solid lines. The northernmost fluorescence detector is called 
Middle Drum while the southern fluorescence detectors are 
referred to as Black Rock Mesa and Long Ridge. The filled 
circle in the middle equally spaced from the three fluorescence 
detectors is the Central Laser Facility used for atmospheric 
monitoring and detector calibration. 


inelastic cross section and on the inclusive proton and 
pion cross sections with the light nuclear atmospheric 
target Q- 

In order to determine AT to derive the interac¬ 
tion length Ap_air from the slope of Xmax distri¬ 
bution Km, we carried out simulation studies using 
the one-dimensional air shower Monte Carlo program 
CONEX4.37( [m, [13, El). The CONEX program uses 
a hybrid air shower calculation for the high energy part of 
the shower, and a numerical solution of the cascade equa¬ 
tions for the low energy part of the shower. This hybrid 
approach of simulating the cosmic ray showers enables 
CONEX to be very efficient. Using CONEX allows us 
to simulate large number of showers in a very reasonable 
time scale. It is worth noting that the shower parameters 
obtained with CONEX are consistent with that obtained 
with CORISKA El- 

Using CONEX the value of AT is determined by sim¬ 
ulating 10,000 events for each of several energy bins for 
data between 10^® ® and 10^® ® eV. The value of K is cal¬ 
culated for each high energy model for each energy bin 
by obtaining the values of Km and Ap_air for that model. 
The value of Km and therefore AT for each of the data 
sets is impacted by the choice of the lower edge of the fit 
range Xi. This dependence is shown in Figure [21 
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It is essential that a consistent procedure be used to 
determine the Xi and consequently the value of K for 
the shower simulations and the observed data. We find 
from the data that Xi =< Xmax > +40 g/cm^ is the 
minimum stable value of Xi, maximizing the number of 
events in the tail of the distribution and consequently the 
statistical power of the measurement. The same relative 
shift distribution is later used in the simulations. 

It is also important to note that in addition to CONEX 
we have also used CORSIKA [l^. CORSIKA is used 
here to simulate three-dimensional cosmic ray showers. 
In the simulation process these showers are thinned in 
order to reduce the CPU time, and then dethinned in an 
attempt to restore lost information M- These show¬ 
ers are then propagated through the FD and the SD 
part of the TA detector. The showers that successfully 
pass the trigger of the detector are then reconstructed, 
after which the pattern recognition event selection are 
applied. The value of Am is then determined and, as 
shown in Figure [31 is found to be consistent with that ob¬ 
tained with CONEX (shower simulation not propagated 
through the detector) particularly around the selected 
choice of Xi =< Xmax > +40 g/cm^. This effect will 
also be discussed in Section CYl 

The value of K is calculated for each simulated data 
set between the energies of 10^®'^ and 10^®®. Figured] 
shows K vs. Logio(E(eV)). Note that we have chosen to 
display QGSJETII.4 as an example. The value of K is 
then established by fitting the points from Figure dl to 
a constant. Table U summarizes the high energy models 
used, the value of K obtained for these models. It is also 
worth mentioning that the K value was also calculated 
with QGSJETII.3 and was obtained from this model to 
be consistent with that determined from QGSJETII.4 
within the statistical fluctuations. Note that the stability 
of K around the average shown in Figuredlshows that K 
is independent of energy and justifies the use of a single 
average value over the range of interest. 

To confirm that the value of K obtained is valid to 
reproduce the interaction length of the model, a plot of 
Ap_air vs. Logio(E(eV)) is shown in Figure|5l Each point 
here represents 10,000 simulated data sets at that energy. 
The circle markers are the Ap_air obtained from the Xi 
distributions from the model, while the triangle markers 
are the Ap_air obtained from reconstruction using the K- 
Factor method. Figuredlshows that using the AT-Factor 
method does indeed reconstruct the expected values of 
the Ap_air for these simulations. This ensures that the 
value of K obtained in this study describe the value of 
K of the high energy models correctly. 

The AT-Factor determined in the procedure described 
above is dependent on the hadronic interaction model 
used in the air shower Monte Garlo simulation. The high 
energy models used in this study are QGSJETII.4 [l6l| . 
QGSJETOl [13, SIBYLL (H, and EPOS-LHC [I^. The 
resultant values of K determined for these models are 
summarized in Table HI 



FIG. 2. The value of K vs. the lower edge in the fit range 
Xi to the tail of the Xmax distribution for several data sets 
10^®'^, 10^®’^, and 10^® eV simulated using CONEX with the 
high energy model QGSJETII.4. Each data set contains 
10,000 simulated events. 
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FIG. 3. Am (g/cm^) vs. the lower edge in the fit range 
Xi to the tail of the Xmax distribution at an energy range of 
10^®'® —eV. The value of Am is calculated using GONEX 
with the high energy model QGSJETII.4 (square markers). 
These events were not propagated through the detector. In 
addition, the value of is also calculated using CORSIKA 
(circle markers). These events successfully survived the pat¬ 
tern recognition cuts after they were successfully detected and 
reconstructed. 


The hrst measurement of the proton-air cross section 
using UHECR was performed by the Flys Eye experi¬ 
ment, which used a calculated value of A" = 1.6 and 
obtained CTp’lfair = 530 ± 66 mb (3 • Following the Fly’s 
Eye result, the calculated values of AT which appeared in 
the literature showed a continuous decrease as full Monte 
Garlo simulations came into use. By 2000, after the de¬ 
velopment of modern high energy hadronic models, the 
reported AT-values still differed by approximately 7% [^ . 
Since then, as shown in Table IH more complete hadronic 
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FIG. 4. The value of K obtained vs. energy in Logic (eV) 
for simulated data sets using CONEX with the high energy 
model QGSJETII.4, for the energy range of the data, between 
10^®-® and 10^®-® eV. 
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FIG. 5. The proton-air interaction length Ap_air in g/cm® 
vs. Energy in Logic (eV) for the simulated data sets using 
GONEX with the high energy model QGSJETIL4, for the 
energy range of the data, between 10^®'® and 10^®'® eV. The 
circle points are the Ap-air values obtained from the Xi dis¬ 
tribution. Triangle points are the ones determined from re¬ 
constructing the Ap_air values using the if-Eactor method. 


Model 

K 

QGSJETIL4 

1.15± 0.01 

QGSJETOl 

1.22T0.01 

SIBYLL 

l.lSiO.Ol 

EPOS-LHG 

1.19T0.01 


TABLE 1. The value of K obtained for each of the high energy 
models. Each K listed is the single average value of K over 
the energy range of 10®® ®-10®® ®. Note that the values of K 
shows a ~3% model uncertainty. 


shower simulations have converged on a smaller value of 
K = 1.2, with a model uncertainty of approximately 3%. 
Using this lower K-value, the Flys Eye cross section may 
be updated to 392 ± 49 mb. 


Model 

o-p^!air [Stat.] mb 

QGSJETIL4 

550.3±68.5 

QGSJETOl 

583.7±72.6 

SIBYLL 

564.6T70.2 

EPOS-LHG 

569.4±70.8 


TABLE 11. The high energy model vs. the CpULir in mb 
obtained for that high energy model. 


IV. PROTON-AIR CROSS SECTION: 

The data used in this analysis is the Telescope Array 
Middle Drum-Surface Detector hybrid events discussed 
in detail in Section |TT1 Figure [6] shows the Xmax distribu¬ 
tion together with the exponential unbinned maximum 
likelihood fit to the tail between 790 and 1000 g/cm^, 
the Am value from the fit is found to be (50.47 ± 
6.26[Stat.]) g/cm^. 

Consecutively the value of CTp^Cj. is determined where 
CTp^fCr = iF X 24,160/Am mb using Equation [TJ The K 
values used are the ones calculated and summarized in 
the previous section in Table [H Accordingly the values of 
^p^air for all the considered hadronic interaction models 
are determined and tabulated in Table |TT1 The final value 
of the proton-air cross inelastic section reported by the 
Telescope Array collaboration is the average value of the 
'^p^air obtained by the high energy models QGJSETII.4, 
QGSJETOl, SIBYLL, and EPOS-LHG and is found to 
be equal to (567.0 ± 70.5[Stat.]) mb. 



FIG. 6. The number of events per Xmax bin (AXmax) vs. 
Xmax g/cm® for the Telescope Array data with the energy 
between 10®®’® and 10®®’® eV. The line is the exponential fit 
to the slope. 

In order to quantify the systematic uncertainties on 
the proton-air cross section obtained using the iF-Factor 
method a few different checks were applied. First, the 
systematic value from the hadronic interaction model 
dependence of the cTp^hj. value is calculated to be the 
maximum difference between the crj/fh,. value determined 
from the various tested models and the average value ob¬ 
tained from these models. The systematic uncertainty 
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from the model dependence is found to be (±17) mb. 

In addition, the systematic error in (Tp“bj. from the 
systematic error in is also calculated. The data is 
divided in halves based on the zenith angle of the events, 
the distance of the shower using the impact parame¬ 
ter, and finally the energy of the events. The attenu¬ 
ation lengths resulting from all these subsets are consis¬ 
tent within the statistical fluctuations. In the case of 
the energy dependence, below the median of 10^® ®^ eV 
A„i = 55.7± 10.1, and above the median A^ = 45.5±7.7. 

Moreover, the systematic effect of possible energy de¬ 
pendent bias in the X^ax distribution was studied. This 
is done by shifting the values of Xmax by their elonga¬ 
tion rate prior to fitting. The value of Am is calculated 
and the systematic effect from a possible energy bias was 
found to be negligible. 

The next check is calculating the systematic uncer¬ 
tainty that originates from the detector bias. This in¬ 
cludes the bias that occurs from detecting the events, 
reconstructing the events, and applying the needed cuts 
to the events. This check is investigated by comparing 
the result of the attenuation length Am of the simulated 
shower thrown without any detector effects to the atten¬ 
uation length obtained from a three-dimensional shower 
simulation using CORSIKA propagated through the de¬ 
tector and reconstructed successfully including the pat¬ 
tern recognition cuts. As shown in Figure [31 the value 
of Am was found to be consistent, for all the high energy 
models, between the thrown events and the reconstructed 
events with pattern recognition applied. Therefore, the 
detector bias systematic effect on the Am value is negli¬ 
gible. 

In addition, a fraction of the high energy cosmic rays 
detected and used in this study are possibly photons. 
Such photons may accompany the cosmic rays by some 
scenarios explaining the origin of ultra high energy cos¬ 
mic ray sources. In addition, a flux of photons is also ex¬ 
pected from the interaction of cosmic rays with energies 
above 4 x lO^^eV with the microwave background radia¬ 
tion produ cing the Greisen-Zatsepin-Kuzmin (GZK) pro¬ 
cess ( [23 )■ There have been several studies placing 

an upper limit on the integral flux and the fraction of 
the primary cosmic ray photons for energies greater than 
10^®eV ( [^, [l^, dll)' this study, the lowest derived 
limit on the photon fraction is used and is < 1% [^ . 
The systematic contribution from the photons is found 
to be -1-23 mb. 

The result of the proton-air cross section from this 
work so far assumes with high energy model simulations 
a pure protonic cosmic ray composition. Regardless of 
what conclusion one would make on the composition in 
the data( the result on the proton-air cross sec¬ 

tion from this work would remain the same. However, the 
systematic effect of the presence of other elements in the 
data beside proton is also studied. This includes iron, 
helium, and CNO. Note that the maximum systematic 
contribution from these elements was found to be from 
helium (deepest Xmax distribution). Hence, It is the con¬ 


Systematic source 

Systematics (mb) 

Model Dependence 

(±17) 

10% Helium 

- 9 

20% Helium 

- 18 

50% Helium 

- 42 

Gamma 

-b 23 

Summary 
(20% Helium) 

(-25,-b29) 


TABLE III. The systematic source vs. The systematic values 
of that source. 


tribution that is reported in this study. A contribution 
of 10%, 20%, and 50% from helium and the systematic 
error associated with such contribution is reported. For 
a 10% contribution the systematic effect is calculated to 
be —9 mb. Meanwhile, for a 20% and 50% contribu¬ 
tion the systematic effect is determined to be —18 mb 
and —42 mb respectively. The final systematic value, 
conservatively assuming a 20% Helium contamination, is 
calculated by adding in quadrature the systematic values 
and is found to be (-25,-1-29) mb. Table ITTTl summarizes 
the systematic checks for the proton-air cross section, in¬ 
cluding the final systematic value. 

We summarize the result of the our proton-air cross 
section obtained using the A"-Factor method described 
previously together with the systematic checks obtained 
to be equal to 


^p™air = 567.0 ± 70.5[Stat.]l25[Sys.]mb. 


( 2 ) 


This is obtained at an average energy of 10^® ®®eV. The 
result of the proton-air cross section is then compared to 
the results obtained from various emerimental results 

( [H, i, [13, in, [13, [13, i, [13) Figures III 


In addition, the experimental results of the high en¬ 
ergy models (QGSJETH.4, QGSJETOl, SIBYLL, EPOS- 
LHG) cross section predictions are also included. This 
includes the statistical (outer/thinner error bar) and the 
systematic (inner/thicker error bar). 


V. PROTON-PROTON CROSS SECTION 

Prom the TA proton-air cross section result we can de¬ 
termine the total proton-proton cross section. The pro¬ 
cess of inferring from is described in details 

in [13, and [13 • 

The Gp-p is calculated from the measured cross section, 
also known as the inelastic cross section 0 ™%^, using 
both Glauber Eormalism [13 and the relation: 

inel _ total _ el _ ^“3^ 

^p—air ^p—air ^p—air ^p—air 

Where is the total cross section, cr®'_air is i-Fe 

elastic cross section and is the quasi elastic cross 

section. The quasi-elastic cross section corresponds to 
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FIG. 7. The proton-air cross section result of this 

work, including the statistical (outer/thinner) and sys¬ 
tematic (inner/thicker) error bar. The result of this 
work is shown in comparison to other experimental results 

( [ill, 01. [H, [al, [al, [al, 0, H). in addition, 

the high energy models (QGSJETII.4, QGSJETOl, SIBYLL, 
EPOS-LHC) cross section predictions are also shown by solid 
line, fine dashed line, dotted line, and dashed line consecu¬ 
tively. 


scattering processes in which nuclear excitation occurs 
without particle production. 

The relation between the and the Up-p is highly 

dependent on the forward scattering elastic slope B. 



In 


C^O'p-p 

dt 




(4) 


This is shown in the B, ( 7 *'^^' plane in Figure [H Here 
the solid and dotted curves represent a constant value of 
f^p^air reflects the Telescope Array measured value 
and the statistical fluctuations. 

There have been many theories predicting the relation¬ 
ship between B and Up-p. However many of these mod¬ 
els either failed to describe the elastic scattering data, or 
the elastic slop e energy dependence from the Tevatron 
( [13, [ 111 , (Hi). A more updated theory using the single 
pomeron exchange model while describing the Tevatron 
data correctly is not consistent with the Unitarity con¬ 
straint ( 111 , 113 ). Here the unitarity constraint is shown 
by solid grey shaded area in Figure [51 A more recent pre- 
diction is the Block, Halzen, and Stanev (BHS) fit H. It 
is consistent with unitarity while using a QCD inspired 
fit to the pp and pp data from the Tevatron. The dashed 
line in Figure [5] shows the BHS prediction. Here 
is converted to using the BHS fit. The statistical 

and systematic errors in are propagated from the 

^p^air statistical and systematic error calculation. The 
CTp'^p is found to be i 70 l 44 IStat.lt 17 [Sys.] mb. 

The calculated in this work is shown in Fig¬ 

ures 0 compared to previous results from cosmic ray ex- 



FIG. 8 . The elastic slope B in ((GeV/c)“^) vs. CTp'llt in 
mb. The solid and the dotted curves are the relation between 
B and for the constant value of the measured 

by the Telescope Array detector and the statistical error us¬ 
ing Glauber Formalism. The dashed line is the BHS QGD 
inspired fit 0. While the gray shaded area is the unitarity 
constraint. 


periments like Fly’s Eye 0, Akeno [13, HiRes [H, and 
Auger [1 , together with accelerator pp and pp cross sec¬ 
tion measurement [dlj, in addition to the recent result 
from LHC by TOTEM [il. The dotted curve is the 
QCD inspired fit of the total p-p cross section vs. the 
center of mass energy ^^(GeV) 0. The result from this 
work at ^/s = 95lg TeV, in addition to the most recent 
result published by the Auger experiment at ^/s = 57 
TeV [I and reported recent result by the LHC at y/s = 
7 TeV [11 are all in agreement with the fit. 
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FIG. 9. The proton-proton cross section vs. the center 
of mass energy result of this work, including the statisti¬ 
cal (outer/thinner) and systematic (inner/thicker) error bars. 
The pp and the pp data are shown in smaller darker circles 
and square symbols consecutively [11 . The recent result from 
LHG is also shown by the star marker [H. The result of this 
work is shown in comparison to previous work by cosmic rays 
detectors( 0, [13, ^], 0). The dashed curve is the QCD 
inspired fit by BHS U 










VI. CONCLUSION AND OUTLOOK 


In this work we used events collected by Telescope Ar¬ 
ray between May-2008 and May-2013 in hybrid mode to 
determine the CTp^aij. using the AT-Factor method. The 
hadronic model dependence of the AT-Factor method was 
investigated. The latest updated hadronic interaction 
models have converged with time on the value of AT 
with an uncertainty of ~3%. This makes the AT-Factor 
method a weakly model dependent method to use in cal¬ 
culating the CTp'lfbj.. Several systematic checks were ap¬ 
plied and the final value of o'p'lfbj, was found to be equal 
to 567.0 ± 70.5[Stat.]t25[Sys.] mb. 

Ultimately the value of crp_p is determined from CTp^fb^ 
using Glauber theory and BHS QCD inspired fit. Such 
a fundamental measurement at this high energy (-^5 = 
95lg TeV) could not be obtained with current parti¬ 
cle accelerators.The value of o'p'^p was determined to be 
170t44[Stat.]l}7[Sys.] mb. 

While the events used in this analysis were collected 
with MD-SD part of the detector, future cross section 
results, using thoroughly analysed events could be per¬ 
formed using LR and BRM data, and ultimately with the 
full detector. LR and BRM are the fluorescence detectors 
closer in distance to the SD and therefore we could extend 
the energy range of the collected data down to 1 EeV . 
This will enable us to study the measurement down to 
1 EeV with higher statistical power which would allow us 
to constrain the available high energy model cross section 
predictions. 
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